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ABSTRACT 

We have analysed the parameters characterising the mass and velocity distributions 
f*^) i of two samples of relaxed elliptical-like-objects (ELOs) identified, at z = 0, in a set of 

£f) ' self-consistent hydrodynamical simulations operating in the context of a concordance 

cosmological model. ELOs have a prominent, non-rotating, dynamically relaxed stellar 
C*~) ■ spheroidal component, with very low cold gas content, and sizes of no more than ~ 10 

- 40 kpc (ELO or baryonic object scale), embedded in a massive halo of dark matter 
typically ten times larger in size (halo scale). They have also an extended halo of 
hot diffuse gas. The parameters characterising the mass, size and velocity dispersion 
both at the baryonic object and at the halo scales have been measured in the ELOs 
of each sample. At the halo scale they have been found to satisfy virial relations; at 
' the scale of the baryonic object the (logarithms of the) ELO stellar masses, projected 

, stellar half-mass radii, and stellar central l.o.s. velocity dispersions define a flattened 

ellipsoid close to a plane (the intrinsic dynamical plane, IDP), tilted relative to the 
virial one, whose observational manifestation is the observed FP. Otherwise, IDPs 
are not homogeneously populated, but ELOs, as well as elliptical (E) galaxies in the 
FP, occupy only a particular region defined by the range of their masses. The ELO 
^—j | samples have been found to show systematic trends with the mass scale in both, the 

relative content and the relative distributions of the baryonic and the dark mass ELO 
components, so that homology is broken in the spatial mass distribution (resulting in 
the IDP tilt), but ELOs are still a two-parameter family where the two parameters 
are correlated (causing its non-homogeneous population) . The physical origin of these 
' trends presumably lies in the systematic decrease, with increasing ELO mass, of the 

relative amount of dissipation experienced by the baryonic mass component along 
ELO stellar mass assembly. ELOs also show kinematical segregation, but it does not 
appreciably change with the mass scale. 

The non-homogeneous population of IDPs explains the role played by M v - lr to 
determine the correlations among intrinsic parameters. In this paper we also show that 
the central stellar line-of-sight velocity dispersion of ELOs, crf*g r , is a fair empirical 
estimator of M v - lri and this explains the central role played by cti OSi o at determining 
the observational correlations. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: haloes - galaxies: kine- 
matics and dynamics - galaxies: structure - dark matter - hydrodynamics 
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1 INTRODUCTION 

Understanding how local galaxies of different Hubble types 
we observe to-day have formed is one of the most chal- 

* Current address: Dept. of Physics, Mahidol University, 
Bangkok 10400, Thailand 



lenging open problems in cosmology. Among the different 
galaxy families, elliptical (Es) are the easiest to study and 
those that show the most precise regularities in their empir- 
ical properties, some times in the form of tight correlations 
among their observable parameters. The interest of these 
regularities lies in that they could encode a lot of relevant 
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informations on the physical processes underlying E forma- 
tion and evolution. 

The Sloan digital sky survey (SDSS, York et al. 2000) 
has substantially improved the statistics on E samples. The 
sample selected by Bernardi et al. (2003a), using morpho- 
logical and spectral criteria, contains 9000 Es to date in the 
redshift range 0.01 < z < 0.3 and in every environment from 
voids to groups to rich clusters. Analyses of their structural 
and dynamical parameters have shown that the distributions 
of their luminosities L, radii at half projected light i? e Ig *, 
and central line-of-sight velocity dispersions <7i os ,o (Bernardi 
et al. 2003b, 2003c), are approximately gaussian at any z. 
Moreover, a maximum likelihood analysis indicates that the 
pairs of parameters o"i os ,o — L and 7?e ght — L, or their com- 
binations, such as the mass-to-luminosity ratio within the 
effective radii M c /L and L (where M e is the dynamical mass 
defined as M c = 2i?' lsht of os /G), show correlations consis- 
tent with those previously established in literature, obtained 
from individual galaxy spectra of smaller samples, such as 
the Faber- Jackson relation (1976); the D n — <7i O s,0 relation 
(Dressier et al. 1987); and the surface brightness — 7? e rela- 
tion (Kormendy 1977, Kormendy & Djorjovski 1989), among 
others. Furthermore, early-type galaxies in the SDSS have 
been found to have roughly constant stellar-mass-to-light 
ratios (Kauffmann et al. 2003a, 2003b; Padmanabhan et al. 
2004). 

The correlations involving two variables out of the three 
needed to fully characterise the structure and dynamics of 
an E galaxy, are projections of the so-called Fundamental 
Plane relation (FP, Djorgovski & Davis 1987; Dressier et 
al. 1987a; Faber et al. 1987; Kormendy & Djorgovski 1989), 
involving the three or some or their combinations. The FP 
relation can be written as 

k>gio R e Sht = ol °gio °"i°s,o + &logio < /hght ><= + c - (!) 

where < J llght > c is the mean surface brightness within 
i^ ight . The values of the FP coefficients for the SDSS E sam- 
ple are a ~ 1.5, similar in the four SDSS bands, b ~ —0.77, 
and c ~ —8.7 (see their exact values in Bernardi et al. 2003c, 
Table 2) with a small scatter. These SDSS results confirm 
previous ones, either in the optical (Lucey, Bower & El- 
lis 1991; de Carvalho & Djorgovski 1992; Bender, Burstein 
& Faber 1992; J0rgensen et al. 1993; Prugniel & Simien 
1996; J0rgensen et al. 1996) or in the near-IR wavelengths 
(Recillas-Cruz et al. 1990, 1991; Pahre, Djorgovski & de Car- 
valho 1995; Mobasher et al. 1999), even if the published val- 
ues of a show larger values in the if-band than at shorter 
wavelengths (see, for example, Pahre, de Carvalho & Djor- 
govski 1998). 

The existence of the FP and its small scatter has the 
important implication that it provides us with a strong con- 
straint when studying elliptical galaxy formation and evo- 
lution (Bender, Burstein & Faber 1993; Guzman, Lucey & 
Bower 1993; Renzini & Ciotti 1993). The physical origin 
of the FP is not yet clear, but it must be a consequence 
of the physical processes responsible for galaxy assembly. 
These processes built up early type galaxies as dynamically 
hot systems whose configuration in phase space are close 
to equilibrium. Taking an elliptical galaxy as a system in 
equilibrium, the virial theorem 



where M v i r is its virial mass, is the average 3- 

dimensional velocity dispersion of the whole elliptical, in- 
cluding both dark and baryonic matter, the dynamical 
half-radius or radius enclosing half the total mass of the 
system, and cf a form factor of order unity, would imply a 
relation similar to Eq. 1 with a = 2 and b = — 1 (known 
as the virial plane, Faber et al. 1987), should the dynamical 
mass-to- light ratios, M v i r /L, and the mass structure coeffi- 
cients 
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M vir = c F (a^) 2 r^/G, 
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be independent of the E luminosity or mass scale. The 
observational results described above mean that the FP is 
tilted relative to the virial plane, and, consequently that ei- 
ther M v i T /L or , or both, do depend on the E luminosity. 

Different authors interpret the tilt of the FP relative 
to the virial relation as caused by different misassumptions 
that we comment briefly (note that we can write M v i T /L = 
M star /L x M vir /M star , where M star is the stellar mass of the 
elliptical galaxy). 

1.1) A first possibility is that the tilt is due to system- 
atic changes of stellar age and metallicity with galaxy mass, 
or, even, to changes of the slope of the stellar initial mass 
function (hereafter, IMF) with galaxy mass, resulting in sys- 
tematic changes in the sieZ/ar-mass-to-light ratios, M star /L, 
with mass or luminosity (Zepf & Silk 1996; Pahre at al. 
1998; Mobasher et al. 1999). But these effects could explain 
at most only ~ one third of the (3 / value in the B- 
band (Tinsley 1978; Dressier et al. 1987; Prugniel & Simien 
1996; see also Renzini & Ciotti 1993; Trujillo, Burkert & 
Bell 2004). Furthermore, early- type galaxies in the SDSS 
have been found to have roughly constant stellar-mass-to- 
light ratios (Kauffmann et al. 2003a, 2003b). Anyhow, the 
presence of a tilt in the JC-band FP, where population effects 
are no important, indicates that it is very difficult that the 
tilt is caused by stellar physics processes alone, as Bender et 
al. (1992), Renzini & Ciotti (1993), Guzman et al. (1993), 
Pahre et al. (1998), among other authors, have suggested. 

1.2) A second possibility is that M v i r /L changes system- 
atically with the mass scale because the total dark-to- visible 
mass ratio, M v i T /M staT changes (see, for example, Renzini 
& Ciotti 1993; Pahre et al. 1998; Ciotti, Lanzoni & Renzini 
1996; Padmanabhan et al. 2004). 

Otherwise, a dependence of cj^ r on the mass scale could 
be caused by systematic differences in: 

2.1) the dark versus bright matter spatial distribution, 

2.2) the kinematical segregation, the rotational support 
and/or velocity dispersion anisotropy in the stellar compo- 
nent (dynamical non- homology) , 

2.3) systematic projection or other geometrical effects. 
Taking into account these effects in the FP tilt demands 

modelling the galaxy mass and velocity three-dimensional 
distributions and comparing the outputs with high quality 
data. 

Bender et al. (1992) considered effects iii) and iv); Ciotti 
et al. (1996) explore ii) - iv) and conclude that an system- 
atic increase in the dark matter content with mass, or dif- 
ferences in its distribution, as well as a dependence of the 
Sersic (1968) shape parameter for the luminosity profiles 
with mass, may by themselves formally produce the tilt; 
Padmanabhan et al. (2004) find evidence of effect ii) in SDSS 
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data. Other authors have also shown that allowing for bro- 
ken homology, either dynamical (Busarello et al. 1997), in 
the luminosity profiles (Trujillo et al. 2004), or both (Prug- 
niel & Simien 1997; Graham & Colless 1997; Pahre et al. 
1998), brings the observed FP closer to Eq. (3). 

Observational methods suffer from some drawbacks to 
study in deph the physical origin of the FP tilt. For exam- 
ple, a drawback is the impossibility to get accurate measure- 
ments of the elliptical three-dimensional mass distributions 
(either dark, stellar or gaseous) and another is that the in- 
trinsic 3D velocity distribution of galaxies is severely limited 
by projection. Only the line-of-sight velocity distributions 
can be inferred from galaxy spectra. And, so, the interpreta- 
tion of observational data is not always straightforward. To 
complement the informations provided by data and circum- 
vent these drawbacks, analytical modelling is largely used in 
literature (Kronawitter et al. 2000; Gerhard et al. 2001; Ro- 
manowsky & Kochanek 2001; Borriello et al. 2003; Padman- 
abhan et al. 2004; Mamon & Lockas 2005a, 2005b). These 
methods give very interesting insights into mass and veloc- 
ity distributions, as well as the physical processes causing 
them, but are somewhat limited by symmetry considerations 
and other necessary simplifying hypotheses. These difficul- 
ties and limitations could be circumvented should we have 
at our disposal complete informations on the phase-space 
of the galaxy constituents. This is not possible through ob- 
servations, but can be attained, at least in a virtual sense, 
through numerical simulations. 

Capelato, de Carvalho & Carlberg (1995) first ad- 
dressed the origin of the FP through numerical simulations. 
By analyzing the remnants of the dissipationless mergers 
of two equal-mass one-component King models, and vary- 
ing their relative orbital energy and angular momentum, 
they show that the mergers of objects in the FP produces a 
new objects in the FP. This result was extended by Dantas 
et al. (2003), who used one- and two-component Hernquist 
models as progenitors, Gozalez- Garcia & van Albada (2003), 
based on Jaffe (1983) models and by Boylan-Kolchin et al. 
(2005), who used Hernquist+NFW models. Nipoti, Londrillo 
& Ciotti (2003) show, in turn, that the FP is well reproduced 
by dissipationless hierarchical equal-mass merging of one- 
and two-component galaxy models, and by accretion with 
substantial angular momentum, with the merging zeroth- 
order generation placed at the FP itself. They also found 
that both the Faber-Jackson and the Kormendy relations 
are not reproduced by the simulations, and conclude that 
dissipation must be a basic ingredient in elliptical formation. 
In agreement with this conclusion, Dantas et al. (2002) and 
Dantas et al. (2003) have shown that the end products of 
dissipationless collapse generally do not follow a FP-like cor- 
relation. Bekki (1998) first considered the role of dissipation 
in elliptical formation through pre-prepared simulations. He 
adopts the merger hypothesis (i.e.., ellipticals form by the 
mergers of two equal-mass gas-rich spirals) and he focuses 
on the role of the timescale for star formation in determin- 
ing the structural and kinematical properties of the merger 
remnants. He concludes that the slope of the FP reflects 
the difference in the amount of dissipation the merger end 
products have experienced according with their luminosity 
(or mass). Recently, Robertson et al. (2006) have confirmed 
this conclusion on the role of dissipative dynamics to shape 
the FP, again through pre-prepared mergers of disk galaxies. 



In this paper we go a step further and we analyse the 
tilt of the FP in samples of virtual ellipticals formed in a 
cosmological context through self-consistent hydrodynami- 
cal simulations. This numerical approach provide a conve- 
nient method to work out the clues of regularity and system- 
atics of elliptical galaxies, and to find out their links with 
the processes involved in galaxy assembly in a cosmologi- 
cal context. The point important for our present purposes 
is that they directly provide with complete 6-dimensional 
phase-space informations on each constituent particle sam- 
pling a given galaxy-like object formed in the simulation, 
that is, they give directly the mass and velocity distribu- 
tions of dark matter, gas and stars of each object. 

Taking advantage of these possibilities, we have anal- 
ysed ten self-consistent hydrodynamical simulations run in 
the framework of a flat ACDM cosmological model, charac- 
terised by cosmological parameters consistent with their last 
determinations (Spergel et al. 2006). Galaxy-like objects of 
different morphologies form in these simulations at z = 0: 
irregulars, disc-like objects, SO-like objects and elliptical- 
like objects (hereafter, ELOs). ELOs have been identified 
as those objects having a prominent dynamically relaxed, 
roughly non-rotating stellar spheroidal component, with no 
extended disks and very low cold gas content; the stellar 
component has typical sizes of no more than ~ 10 - 40 
kpc, it dominates the mass density at these scales (here- 
after, ELO or baryonic object scale) and it is embedded in a 
halo of dark matter typically ten times larger in size (here- 
after, halo scale). In a forthcoming paper (Onorbe et al., in 
preparation) we report on an analysis of the mass and ve- 
locity distributions of the different ELO components: dark 
matter, stars, cold gas and hot gas. In this paper we focus on 
the quantitative characterisation of these mass and velocity 
distributions through their corresponding parameters, both 
at the ELO and at the halo scales. At the baryonic object 
scale, to characterise the structural and dynamical proper- 
ties of ELOs, we will describe their three dimensional dis- 
tributions of mass and velocity through the three intrinsic 
(that is, three-dimensional) parameters (the stellar mass at 
the baryonic object scale, M^ T , the stellar half- mass radii 
at the baryonic object scale, rjfbo, defined as those radii en- 
closing half the M£o mass, and the mean square velocity 
for stars, o-f ^ T ) whose observational projected counterparts 
(the luminosity L, effective projected size i^ ight , and the 
stellar central l.o.s. velocity dispersion, cti os ,o) enter the def- 
inition of the observed FP (Eq. 1). To help the reader, in 
Table 1 we give a list of the parameter names and symbols 
and in Table 2 a list of the profiles and ratios. These in- 
trinsic parameters have been measured on ELOs, and their 
correlations have been looked for, and more specifically, the 
lower-dimensional regions (i.e., dynamical planes) they fill in 
the three dimensional space of the M{f ar , rf£o and °1i!bo Pa- 
rameters, because the observational manifestation of these 
dynamical planes is the FP relation. These informations, 
combined with that on the mass, size and velocity disper- 
sion parameters at the halo scale, allows us to test whether 
or not the Cm* coefficients and the M v i r /M star ratios do sys- 
tematically depend on the mass scale, so that the tilt and 
the scatter of the observed FP can be explained in terms of 
the regularities in the structural and dynamical properties of 
ELOs formed in self-consistent hydrodynamical simulations. 

Fully-consistent gravo-hydrodynamical simulations as a 
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method to study E assembly has already proven to be use- 
ful. An analysis of ELO structural and kinematical proper- 
ties that can be constrained from observations (i.e., stellar 
masses, projected half-mass radii, central line-of-sight veloc- 
ity dispersions), has shown that they have counterparts in 
the local Universe as far as these properties are concerned 
(Saiz 2003; Saiz, Dommguez-Tenreiro & Serna 2004, here- 
after SDTS04), including the FP relation and some clues 
about its physical origin (Onorbe et al. 2005), and its lack of 
dynamical evolution (Dommguez-Tenreiro et al. 2006, here- 
after DTal06). Also, ELO stellar populations have age distri- 
butions showing similar trends as those inferred from obser- 
vations (Dommguez-Tenreiro, Saiz & Serna 2004, hereafter 
DSS04). 

The paper is organised as follows: in §2 we briefly de- 
scribe the simulations, the ELO samples and their generic 
properties. The ELO size and mass scales and their relations 
are analysed in §3. §4 is devoted to kinematics and in §5 
we report on the intrinsic dynamical plane of ELOs and its 
comparison with the observed Fundamental Plane, and we 
analyse its physical origin. Finally, in §6 we summarise our 
results and discuss them in the context of theoretical results 
on halo structure and dissipation of the gaseous component. 



2 THE SIMULATIONS AND THE ELO 
SAMPLES 

We have analysed ELOs identified in ten self-consistent cos- 
mological simulations run in the framework of the same 
global flat ACDM cosmological model, with h = 0.65, 
Sl m = 0.35, fib = 0.06. The normalisation parameter has 
been taken slightly high, as = 1.18, as compared with the 
average fluctuations of 2dFGRS or SDSS galaxies (Lahav et 
al. 2002, Tegmark et al. 2004) or recent results from WMAP 
(Spergel et al. 2006) to mimic an active region of the Uni- 
verse (Evrard, Silk & Szalay 1990). 

We have used a lagrangian code, (DEVA, Serna, 
Dommguez-Tenreiro & Saiz, 2003) , particularly designed to 
study galaxy assembly in a cosmological context. Gravity is 
computed through an AP3M-like method, based on Couch- 
man (1991). Hydrodynamics is computed through a SPH 
technique where special attention has been paid to make the 
implementation of conservation laws (energy, entropy and 
angular momentum) as accurate as possible (see Serna et al. 
2003 for details, in particular for a discussion on the obser- 
vational implications of violating some conservation laws). 
Entropy conservation is assured by taking into considera- 
tion the space variation of the smoothing length (i.e., the 
so-called V/i terms). Time steps are individual for particles 
(to save CPU time, allowing a good time resolution) , as well 
as masses. Time integration uses a PEC scheme. In any run, 
an homogeneously sampled periodic box of 10 Mpc side has 
been employed and 64 3 dark matter and 64 3 gas particles, 
with a mass of 1.29x 10 8 and 2.67x 10 7 M©, respectively, have 
been used. The gravitational softening used was e = 2.3 kpc. 
The cooling function is that from Tucker (1975) and Bond 
et al. (1984) for an optically thin primordial mixture of H 
and He (X — 0.76, Y = 0.24) in collisional equilibrium and 
in absence of any significant background radiation field with 
a primordial gas composition. Each of the ten simulations 
started at a redshift z m = 20. 



Star formation (SF) processes have been included 
through a simple phenomenological parametrisation, as that 
first used by Katz (1992, see also Tissera et al. 1997; Saiz 
2003 and Serna et al. 2003 for details) that transforms 
cold locally-collapsing gas at the scales the code resolves, 
denser than a threshold density, pthrcs, into stars at a rate 
dpstar/Gfe = c»p gas /fg, where t g is a characteristic time- 
scale chosen to be equal to the maximum of the local gas- 
dynamical time, tdyn = (47rGpgas) -1 ^ 2 , and the local cooling 
time; c* is the average star formation efficiency at resolution 
e scales, i.e., the empirical Kennicutt- Schmidt law (Kenni- 
cutt 1998). It is worth noting that, in the context of the 
new sequential multi-scale SF scenarios (Vazquez-Semadeni 
2004a, 2004b; Ballesteros-Paredes et al. 2006 and references 
therein), it has been argued that this law, and particularly 
so the low c* values inferred from observations, can be ex- 
plained as a result of SF processes acting on dense molecular 
cloud core scales when conveniently averaged on disc scales 
(Elmegreen 2002; Sarson et al. 2004, see below). Supernova 
feedback effects or energy inputs other than gravitational 
have not been explicitly included in these simulations. We 
note that the role of discrete stellar energy sources at the 
scales resolved in this work is not yet clear, as some authors 
argue that stellar energy releases drive the structuring of the 
gas density locally at sub- kpc scales (Elmegreen 2002). In 
fact, some MHD simulations of self-regulating SNII heating 
in the ISM at scales < 250 pc (Sarson et al. 2004), indicate 
that this process produces a Kennicutt-Schmidt-like law on 
average. If this were the case, the Kennicutt-Schmidt law 
implemented in our code would already implicitly account 
for the effects stellar self-regulation has on the scales our 
code resolves, and our ignorance on sub-kpc scale processes 
would be contained in the particular values of pthrcs and c, . 

Five out of the ten simulations (the SF-A type simula- 
tions) share the SF parameters (pthrcs = 6 x 10~ 25 gr cm~ 3 , 
c* = 0.3) and differ in the seed used to build up the initial 
conditions. To test the role of SF parameterisation, the same 
initial conditions have been run with different SF parame- 
ters (pthros = 1-8 x 10~ 24 gr cm~ 3 , c = 0.1) making SF 
more difficult, contributing another set of five simulations 
(hereafter, the SF-B type simulations). 

ELOs have been identified in the simulations as those 
galaxy-like objects at z = having a prominent, non- 
rotating, dynamically relaxed spheroidal component made 
out of stars, with no extended discs and very low cold gas 
content. It turns out that, at z — 0, 26 (17) objects out of the 
more massive formed in SF-A (SF-B) type simulations fulfil 
this identification criterion, forming two samples (the SF-A 
and SF-B ELO samples) partially analysed in SDTS04, in 
DSS04 and in DTal06. In Onorbe et al. (2005) it is shown 
that both samples satisfy dynamical FP relations. ELOs in 
the SF-B sample tend to be of later type than their corre- 
sponding SF-A counterparts because forming stars becomes 
more difficult; this is why many of the SF-B sample coun- 
terparts of the less massive ELOs in SF-A sample do not 
satisfy the selection criteria, and the SF-B sample has a 
lower number of ELOs that the SF-A sample. 

A visual examination of ELOs indicates that the stellar 
component is embedded in a dark matter halo, contribut- 
ing an important fraction of the mass at distances from the 
ELO centre larger than ~ 15 kpc on average. ELOs have 
also a hot, extended, X-ray emitting halo of diffuse gas (Saiz, 
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Dommguez-Tenreiro & Serna 2003). Stellar and dark mat- 
ter particles constitute a dynamically hot component with 
an important velocity dispersion, and, except in the very 
central regions, a positive anisotropy. In Table 3 and Ta- 
ble 4 different data on ELOs in the samples are given. In 
these Tables, the criterion introduced in Saiz et al. 2001, 
ELOs have been labelled by a three-digit code formed by 
rounding their x, y and z coordinates in units where the 
simulated box has, at any redshift, a length of one. The 
number of particles of each species sampling each ELO in 
the sample can be easily determined from these Tables and 
the values of the masses of dark and baryonic particles. The 
spin parameters of the ELO samples have an average value 
of A = 0.033. ELO mass function is consistent with that of 
a small group environment (Cuesta-Bolao & Serna, private 
communication) . 

The simulations show the physical patterns of ELO 
mass assembly and star formation (SF) rate histories (see 
DSS04 and DTal06 for more details). ELOs form out of the 
mass elements that at high z are enclosed by overdense re- 
gions whose mass scale (total mass they enclose) is of the 
order of an E galaxy total (i.e., including its halo) mass. 
Analytical models, as well as N-body simulations indicate 
that two different phases operate along halo mass assembly: 
first, a violent fast one, where the mass aggregation rates are 
high, and then, a slower one, with lower mass aggregation 
rates (Wechsler et al. 2002; Zhao et al. 2003; Salvador-Sole, 
Manrique, & Solanes 2005). Our hydrodynamical simula- 
tions indicate that the fast phase occurs through a multi- 
clump collapse (see Thomas, Greggio & Bender 1999) ensu- 
ing turnaround of the overdense regions, and it is charac- 
terised by the fast head-on fusions experienced by the nodes 
of the cellular structure these regions enclose, resulting in 
strong shocks and high cooling rates of their gaseous com- 
ponent, and, at the same time, in strong and very fast star 
formation bursts that transform most of the available cold 
gas into stars. Consequently, most of the dissipation involved 
in the mass assembly of a given ELO occurs in this vio- 
lent early phase at high z. The slow phase comes after the 
multiclump collapse. In this phase, the halo mass aggrega- 
tion rate is low and the mass increments result from major 
mergers, minor mergers or continuous accretion. Our cosmo- 
logical simulations show that the fusion rates are generally 
low, and that a strong star formation burst and dissipation 
follow a major merger only if enough gas is still available 
after the early violent phase. This is very unlikely in any 
case, and it becomes more and more unlikely as the ELO 
mass increases (see DSS04). And so, these mergers imply 
only a modest amount of energy dissipation or star forma- 
tion, but they play an important role in this slow phase: 
an ~ 50% of ELOs in the sample have experienced a major 
merger event at 2 < z < 0, that result in the increase of the 
ELO mass content, size, and stellar mean square velocity. 
So, our simulations indicate that most of the stars of to-day 
ellipticals formed at high redshifts, while they are assem- 
bled later on (see de Lucia et al. 2006, for similar results 
from a semi-analytic model of galaxy formation grafted to 
the Millennium Simulation). This scenario shares some char- 
acteristics of previously proposed scenarios (see discussion 
and references in DTal06), but it has also significant differ- 
ences, mainly that most stars form out of cold gas that had 
never been shock heated at the halo virial temperature and 



then formed a disk, as the conventional recipe for galaxy for- 
mation propounds (see discussion in Keres et al. 2005 and 
references therein). 



3 SIZE AND MASS SCALES 

3.1 Masses and Sizes at the Scale of the Virial 
Radius 

The virial radius describes the ELO size at the scale of its 
dark matter halo. It roughly encloses those particles that are 
bound into the self-gravitating configuration forming a given 
ELO system (i.e., a dark matter halo plus the main baryonic 
compact object plus the substructures and satellites hosted 
by the dark matter halo). The virial radii have been cal- 
culated using the Bryan & Norman (1998) fitting formula, 
that yields, at z — 0, a value of A ~ 100 for the mean den- 
sity within r v i r in units of the critical density. The mass at 
the scale of r v i r is the virial mass, M v i r , the total mass in- 
side r v ir or halo total mass. The mass scales associated to 
the different constituents considered here are^ : dark matter, 
M^ ark , baryons of any kind, M^ b , cold baryons (that is, cold 
gas particles with T < 3 x 10 4 K and stellar particles), M£ b , 
stars, Mff ar , and hot gas (that is, gaseous particles with 
T > 3 x 10 4 K). A measure of the compactness of the mass 
distribution for the different ELO constituents, at the halo 
scale, is given by their respective half-mass radii, or radii 
enclosing half the mass of these constituents within r v ; r ; for 
example, the overall half-mass radii, r* ^, are the radii of 
the sphere enclosing M v i r /2, the stellar half-mass radii r^ r 
enclose Mjf ar /2 and so on. 

For each of the ELOs in our sample, their virial mass 
and radii are listed in Table 3, where also given are the 
masses within r v j r corresponding to different constituents 
and some relevant half-mass radii. All these mass scales are 
strongly correlated with M V [ T as shown in Figure 1 for Mjf ar . 
Note that the virial masses of ELOs have a lower limit of 
3.7 x 10 11 M Q . 

An important point is the amount of gas infall rel- 
ative to the halo mass scale. As illustrated in Figure 2 
for M£ b /M vir , any of the ratios M£ b /M vil , M£ b /M vir or 
Mif ar /M v j r decreases as M v i r increases, as observationally 
found at smaller scales (see §1). Note that we have in any 
case M^ b /M v j r < Clt,/Q m = 0.171, the average cosmic frac- 
tion, so that there is a lack of baryons within r v i r relative to 
the dark mass content that becomes more important as M V [ T 
increases. Otherwise, heating processes along ELO assembly 
give rise to a hot gas halo around the objects, partially be- 
yond the virial radii. The amount of hot gas mass outside 
the virial radii, normalised to the ELO stellar mass M^ ar 
(see §3.2), increases with the mass scale. It also increases rel- 
ative to the cold gas content at the halo scale. To illustrate 
this point, in Figure 3 we plot the integrated hot gas den- 
sity profile normalised to the mass of cold baryons inside the 
virial radii, M£ b . The mass effect can be clearly appreciated 
in this Figure, where we see the following: 

(i) The mass of hot gas increases monotonically up to 

t Note that wc have used superscripts to mean the different ELO 
constituents, and subscripts to distinguish between halo (h) or 
baryonic object (bo) scales, see Table 1 
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Table 1. Parameter names and symbols 



[H] 



Name 



Symbol 



Observational parameters 



Luminosity 

Half projected light radius 
Central LOS velocity dispersion 
Dynamical mass 

Mean surface brightness within i?j! sht 
Stellar Mass 

Stellar-mass-to-light ratio 



L 

plight 

Clos.O 
M c 

M star 

^star 



Halo scale parameters 



Virial mass 
Virial radius 

Dark mass inside virial radius 
Baryon mass inside virial radius 
Cold baryon mass inside virial radius 
Stellar mass inside virial radius 
Total half-mass radius 
Cold baryon half-mass radius 
Stellar half-mass radius 
Total 3D velocity dispersion 



M vir 

?Vir 
M dark 

M ab 
M£ b 

M star 

„tot 
Vh 
r .cb 
c,h 
„star 
r e,h 
tot 
°3,h 



Baryonic-object scale parameters 



Stellar mass 

Cold baryon mass 

Stellar half-mass radius 

Cold baryon half-mass radius 

Projected stellar half-mass radius 

Mean stellar 3D velocity dispersion 

Central LOS stellar velocity dispersion 

Mean projected stellar mass density within R B 



M star 

M? h 

bo 
^star 
e.bo 
r cb 

e,bo 
ostar 
e,bo 
_star 
°3,bo 
star 
°los,0 
^star ^ 



Table 2. Profiles and ratios 



Profiles Ratios 



Name 


Symbol 


Ratio definition 


Ratio symbol 


Logarithmic slope 


Hot baryon mass profile 

Circular velocity profile 

3D velocity dispersion profile 

Anisotropy profile 

Projected mass density profile 

Line-of-sight velocity profile 

Line-of-sight velocity dispersion profile 


M hb (r) 

Vcir(r) 

°- 3 dM 

AniM 

S(R) 

Vios(K) 
^los(-R) 


G*vi,/Kt) ! r« 

^tot /^star 

c,h' e,bo 
^.star /pstar 

e,bo' e,bo 
(-.tot /_star \2 
V°3,h/ 3,bo/ 

(<Co) 2 /3«%) 2 

GM vir /3« s a ; ) 2 i?i lght = c P c rdCrpCvd c vpc 


c F 

Crd 
Crp 
Cvd 
Cvpc 

„vir 
C M 


0F 

/3 rd 
A P 
Ad 

/3vpc 

Ai 



(a) To specify the constituent, a superindex has been added in the text to the profile symbols 



r ~ 4r v ir, and maybe also beyond this value, but it is difficult 
at these large radii to properly dilucidate whether or not 
a given hot gas mass element belongs to a given ELO or 
to another close one (to alleviate this difficulty, only those 
ELOs not having massive neighbours within radii of 6xr v i r 
have been considered to draw this Figure) . 



(ii) The hot gas mass fraction increases with M vir at 
given r/rvir. This suggests that the cold baryons that mas- 
sive ELOs miss inside r v i r relative to less massive ones ap- 
pear as a diffuse warm component at the outskirts of their 
configurations. 

We now comment on length scales. The overall half- 
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Table 3. Masses, sizes and mean square velocities of ELOs at the halo scale (z = 0) 



Run ELO M vir M^ ark M* b M£ b M= tar r vir r* * rfg 1 



8714 #173 


772.82 


678.31 


94.51 


70.65 


66.69 


527.00 


222.54 


64.80 


51.63 


302.43 


#353 


322.18 


285.71 


36.47 


31.60 


29.61 


394.00 


110.99 


18.62 


16.11 


261.18 


#581 


177.03 


156.26 


20.77 


18.94 


17.59 


322.00 


101.54 


19.05 


16.07 


202.53 


#296 


153.52 


134.38 


19.14 


16.79 


15.56 


308.00 


92.12 


17.09 


13.20 


201.13 


#373 


74.00 


64.23 


9.77 


8.35 


7.66 


241.00 


66.65 


6.70 


5.79 


162.36 


#772 


60.71 


52.16 


8.55 


7.50 


6.89 


226.00 


56.57 


6.58 


5.40 


159.75 


#284 


53.76 


46.03 


7.73 


6.86 


6.15 


217.00 


67.69 


6.34 


4.96 


149.92 



8747 



#288 
#115 
#189 
#915 



285.17 
178.07 
45.15 
38.63 



251.80 
157.31 
38.62 
32.47 



33.37 
20.77 
6.53 
6.17 



29.83 
17.62 
5.77 
5.51 



28.27 
16.14 
5.42 
5.05 



378.00 
323.00 
205.00 
194.00 



112.14 
95.56 
58.45 
47.74 



25.87 
11.28 
4.88 
4.14 



22.65 
8.24 
4.45 
3.59 



241.90 
220.92 
145.62 
141.02 



8741 #011 


273.34 


237.46 


35.87 


30.12 


29.42 


373.00 


102.20 


21.83 


20.37 


246.70 


#017 


141.44 


121.97 


19.47 


16.33 


15.29 


299.00 


89.04 


17.01 


13.93 


197.56 


#930 


135.71 


118.56 


17.15 


13.70 


12.67 


295.00 


89.82 


20.45 


15.95 


209.05 


#945 


107.34 


93.26 


14.08 


12.63 


11.59 


273.00 


78.06 


11.07 


9.00 


184.78 


#097 


108.18 


95.20 


12.97 


11.33 


10.45 


274.00 


73.36 


9.80 


8.45 


185.69 


#907 


55.15 


47.67 


7.48 


6.35 


5.81 


219.00 


50.61 


4.32 


3.77 


160.98 



8742 



#234 
#283 
#254 
#092 



296.59 
160.41 
147.26 
75.03 



260.31 
137.78 
128.57 
64.57 



36.29 
22.63 
18.69 
10.46 



29.86 
16.60 
16.98 
9.14 



28.80 
16.10 
15.37 
8.46 



383.00 
312.00 
303.00 
242.00 



101.79 
87.69 
100.62 
84.18 



15.71 
9.47 
17.43 
9.84 



14.62 258.68 

8.90 235.18 

12.99 190.07 

8.52 151.06 



8743 #238 


327.12 


289.42 


37.70 


33.03 


30.92 


396 


.00 


116.33 


24.41 


20.78 


255.62 


#328 


66.83 


56.77 


10.06 


8.63 


8.36 


233 


.00 


56.61 


5.64 


5.44 


167.70 


#515 


56.33 


48.75 


7.58 


6.79 


6.40 


220 


.00 


56.20 


4.66 


4.34 


158.78 


#437 


52.65 


45.23 


7.42 


6.78 


6.31 


215 


.00 


54.73 


8.00 


6.83 


152.63 


#421 


36.75 


31.44 


5.31 


4.69 


4.29 


191 


.00 


46.77 


4.33 


3.82 


141.13 



8716 


#173 
#253 
#581 
#296 


753.34 
312.64 
170.89 
153.65 


673.07 
281.54 
151.97 
135.26 


80.28 
31.10 
18.92 
18.39 


64.75 
27.42 
17.26 
16.36 


54.89 
24.03 
14.68 
13.34 


523.00 
390.00 
319.00 
308.00 


230.70 
109.32 
100.06 
93.98 


39.85 
8.85 
8.65 
8.39 


24.19 
6.72 
6.07 
5.07 


300.90 
261.01 
207.37 
205.79 


8717 


#317 
#288 
#348 


739.03 
280.20 
222.99 


673.61 
249.36 
197.69 


65.41 
30.85 
25.30 


62.44 
28.17 
23.15 


56.86 
24.27 
19.23 


519.00 
376.00 
348.00 


157.32 
108.00 
111.59 


30.96 
14.90 
14.26 


24.15 
9.54 
8.81 


345.96 
245.54 
222.23 


8721 


#011 

#945 
#097 


271.20 
107.97 
106.32 


237.67 
94.22 
93.99 


33.53 
13.75 
12.33 


28.36 
12.44 
10.54 


26.61 
10.20 
8.99 


372.00 
274.00 
272.00 


105.34 
78.03 
73.94 


9.83 
5.60 
5.14 


7.90 
3.73 
3.83 


248.27 
187.11 
185.74 


8722 


#293 
#234 
#283 
#254 


729.97 
292.30 
157.20 
145.18 


658.54 
259.61 
136.37 
127.78 


71.42 
32.69 
20.83 
17.40 


65.54 
27.44 
15.78 
16.04 


58.46 
25.72 
13.95 
13.13 


516.00 
381.00 
310.00 
302.00 


151.34 
100.25 
86.88 
102.94 


28.27 
7.96 
4.84 
6.47 


20.95 
6.87 
3.62 
4.33 


332.75 
261.70 
237.45 
193.14 


8723 


#647 
#238 
#563 


773.24 
318.02 
271.83 


689.69 
285.59 
240.05 


83.55 
32.44 
31.78 


71.59 
29.60 
27.93 


59.77 
24.81 
24.49 


527.00 
392.00 
372.00 


182.08 
107.38 
110.79 


92.82 
12.58 
24.03 


44.22 
8.06 
15.38 


325.39 
260.00 
240.78 



Masses are given in 10 10 Mq, distances in kpc, velocity dispersion in km s 1 . 



mass radii r*°^ are closely correlated to M V [ T . Concerning 
baryon mass distributions, dissipation in shocks and gas 
cooling play now important roles to determine these mass 
distributions. And so, the rj! " radii depend on how much 
energy was radiated before gaseous particles became dense 
enough to be turned into stars. This, in turn, depends on 
the mass scale, on the one hand, and, in a given mass range, 



on the values of SF parameters, on the other hand. And so, 
more massive ELOs tend to have larger radii and, in 
a given mass range, SF-A sample ELOs tend to have larger 
rjf" radii than their SF-B sample counterparts, because the 
SF implementation in the code demands denser gas to form 
stars in the later than in the former. This effect is more re- 
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M vlr (M e ) 

Figure 1. Masses at the halo scale of stars versus their corre- 
sponding virial masses, filled black symbols: SF-A sample ELOs; 
open red symbols: SF-B sample ELOs 



M vlr (M e ) 

Figure 2. Masses of cold baryons inside the virial radii in units 
of the corresponding virial mass for the ELO sample, filled black 
symbols: SF-A sample ELOs; open red symbols: SF-B sample 
ELOs. 



markable for sizes at the scale of the baryonic object, as we 
shall see in the next subsection. 




Figure 3. The hot gas density profile normalised to the mass 
of cold baryons inside the virial radii for isolated ELOs. Orange 
point-dashed lines: ELOs in the mass interval 1.5 X 10 12 Mq < 
M v i r < 5 X 10 12 Mq; continuous green lines: ELOs in the mass 
interval M vir < 1.5 X 10 12 Mq. 



is M£q, the total amount of cold baryons that have reached 
the central volume of the haloes, forming an ELO. Most of 
these cold baryons have turned into stars, depending on the 
strength of the dynamical activity in the volume surround- 
ing the proto-ELO at high z, and, also, on the values of the 
SF parameters. M^o is the stellar mass. It can be estimated 
from luminosity data through modelling (see Kauffmann et 
al. 2003a, for example"^"). Effective or half-mass radii at the 
baryonic object scale, rf\, and r^^ a , can be defined as those 
radii enclosing half the M{£ or M^ 1 masses, respectively. 
These are the relevant size scales for the intrinsic ELOs, 
but the observationally relevant size scales are the projected 
half- mass radii. They are determined from M cy \{R), the inte- 
grated projected mass density in concentric cylinders of ra- 
dius R for the different constituents. For example, R% ho and 
Rf^l are the projected radii where M^(R) and M c s ^ r (7?) 
are equal to and to M^ 1 /2, respectively. Significant 

parameters at the baryonic object scale are listed in Table 4. 
Note that ELOs have a lower limit in their stellar mass con- 
tent of 3.8 x 10 10 Mq (see Kauffmann et al. 2003b for a 
similar result in SDSS early-type galaxies). 

To illustrate how the halo total mass, M v i r , determines 
the ELO structure at kpc scales, in Figures 4 and 5 we draw 
Mbo ar and rf ^l versus M vir , respectively, for the ELO sam- 
ple. A good correlation is apparent in Figure 4, where it is 
shown that ELO stellar masses are mainly determined by 
the halo mass scale, M v i r , with only a very slight depen- 
dence on the SF parametrisation (SF-A type ELOs have a 
slightly higher stellar content than their SF-B counterparts, 
as expected) . Figure 5 shows also a good correlation between 



3.2 Masses and Sizes at the Scale of the Baryonic 
Object 

Let us now turn to the study of ELOs at the scale of the 
baryonic objects themselves, that is, at scales of some tens 
of kpcs. Physically, the mass parameter at the ELO scale 



t The results of these authors indicate that for SDSS elliptical 
galaxies the stellar-mass-to-light ratio, -y star J can be taken to be 
constant in the range of absolute luminosities M < —21 (see 
Kauffmann et al. 2003b). The values of the logarithm of this ratio 
are log7^ tar ~ 0.53 and log7| tar ~ 0.25, with dispersions erg < 
0.15 and 0.1, in the r and z SDSS bands, respectively. 
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Table 4. Masses, sizes and mean square velocities at the baryonic 10 
object scale, as well as projected radii and central stellar l.o.s. 
velocity dispersions (z = 0) 



Run 


ELO 


bo 


bo 




e, bo 


"Ij.bo 




8714 


#173 


43.12 


42.59 


13.01 


12.72 


351.93 


226.67 




#353 


25.62 


25.17 


8.25 


7.95 


297.11 


192.47 




#581 


13.72 


13.45 


6.75 


6.66 


222.42 


136.18 




#296 


12.21 


11.95 


4.76 


4.69 


225.95 


138.39 




#373 


7.38 


7.23 


3.85 


3.83 


186.95 


116.88 




#772 


5.62 


5.51 


2.33 


2.30 


181.59 


107.28 




#284 


5.70 


5.48 


2.81 


2.74 


169.09 


103.74 


8747 


#288 


20.11 


20.05 


7.84 


7.81 


269.58 


171.50 




#115 


12.21 


12.17 


2.94 


2.96 


259.05 


161.63 




#189 


4.74 


4.68 


2.57 


2.61 


169.93 


104.81 




#915 


4.52 


4.32 


2.23 


2.19 


164.69 


99.61 


8741 


#011 


19.10 


19.10 


5.36 


5.35 


279.61 


174.62 




#017 


10.34 


10.15 


4.15 


4.09 


224.37 


139.54 




#930 


8.95 


8.85 


4.68 


4.70 


214.97 


138.64 




#945 


10.00 


9.89 


4.58 


4.62 


209.68 


128.84 




#097 


9.46 


9.34 


5.00 


5.05 


206.21 


132.19 




#907 


5.63 


5.55 


2.57 


2.60 


187.67 


119.71 


8742 


#234 


27.42 


27.13 


9.33 


9.23 


287.18 


190.47 




#283 


13.51 


13.44 


4.28 


4.30 


246.54 


151.56 




#254 


10.91 


10.81 


4.50 


4.50 


221.72 


134.23 




#092 


6.55 


6.38 


3.14 


3.17 


177.69 


108.46 


8743 


#238 


20.97 


20.88 


6.76 


6.79 


289.43 


181.53 




#328 


7.76 


7.60 


3.40 


3.36 


187.15 


116.49 




#515 


5.74 


5.63 


2.60 


2.62 


185.04 


114.13 




#437 


4.79 


4.70 


2.36 


2.40 


172.71 


105.13 




#421 


3.91 


3.81 


2.11 


2.14 


161.41 


102.09 


8716 


#173 


39.13 


38.07 


7.07 


6.75 


355.08 


237.07 




#253 


22.98 


22.54 


4.40 


4.28 


305.37 


199.34 




#581 


13.03 


12.46 


3.35 


3.18 


243.71 


152.89 




#296 


10.87 


10.43 


2.52 


2.50 


241.58 


140.14 


8717 


#317 


42.28 


41.55 


8.12 


7.77 


387.30 


223.73 




#288 


18.63 


18.25 


3.48 


3.45 


286.59 


169.81 




#348 


15.36 


14.84 


3.73 


3.57 


261.06 


158.31 


8721 


#011 


21.73 


21.37 


3.34 


3.31 


287.47 


177.04 




#945 


9.99 


9.30 


2.55 


2.44 


226.33 


128.90 




#097 


8.62 


8.34 


2.34 


2.36 


213.00 


133.46 


8722 


#293 


43.56 


42.75 


9.40 


9.09 


365.95 


215.68 




#234 


22.40 


22.28 


3.53 


3.57 


313.03 


196.62 




#283 


12.67 


12.51 


2.10 


2.13 


262.81 


163.89 




#254 


11.97 


11.52 


2.64 


2.62 


237.10 


137.68 


8723 


#647 


34.01 


33.43 


5.86 


5.76 


336.69 


210.25 




#238 


20.09 


19.91 


3.78 


3.74 


297.51 


191.05 




#563 


13.67 


13.44 


2.41 


2.45 


274.84 


164.01 



10 12 10 13 
M vlr (M e ) 

Figure 4. Stellar masses at the baryonic object scale scale versus 
halo mass for the ELO sample. Symbols are as in previous Figures 



the length scales for the stellar masses and M v i r , but now 
the sizes depend also on the SF parameters. The physical 
foundations of this behaviour are the same as discussed in 
§3.1. 

The observationally relevant scalelengths are the pro- 
jected half-mass radii ii^bo- Their correlations with their 
intrinsic three dimensional counterparts rj^o are very good, 
as illustrated in Figure 6, where the very low dispersion in 
the plots of the c rp = rf^/Rf^ ratios versus the stellar 
mass M^ T can be appreciated. The results of a fit to a 
power law of the form c rp = A tp (M^ lt ) i3tip § are given in 
Table 5 , where we see that the c rp ratios show a very mild 
mass dependence in the SF-A sample and none in the SF-B 
sample. This result is important because it indicates that 
the observationally available projected radii -R^bo are robust 
estimators of the physically meaningful size scales rj'jjj^. 

We now address the correlations of normalised mass and 
size scales. The increasing behaviour of the M V xt/Mj£ and 
M v i r /M D o ar ratios with increasing mass scale are very inter- 
esting. In particular, the last ratio (Figure 7) follows the 
same trends as the empirical M c /L versus L relation, see 
Bernardi et al. 2003b. The results of a fit to a power law of 
the form M vir /M^ ar = A vit {M^) ^ are given in Table 5, 
where we see that they do not depend on the SF parame- 
terisation. 

To have an idea on how important cold baryon infall 
has been at the baryonic object scale relative to that at the 
halo scale, in Figure 8 the M£ /M^„ ratios are drawn as a 
function of the ELO mass scale. We see that in any case 



Masses are given in 10 10 M , distances in kpc, velocity disper- more than half the mass of cold baryons inside the virial 

sions in km s _1 . radii are concentrated in the central baryonic object, and 

that there is a mass effect in the sense that this fraction 



S We list the different ratio definitions we use in this paper and 
their corresponding logarithmic slopes in Table 2 
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Table 5. Slopes for Linear Fits 









SF-A 




SF-B 








Mi 


0.256 


± 0.035 


0.281 


± 0.048 






/3vir 


0.221 


± 0.083 


0.237 


± 0.158 






Pm 


-0.204 


± 0.116 


-0.247 


± 0.189 






/3f 


0.025 


± 0.048 


0.022 


± 0.081 








0.021 


± 0.041 


0.076 


± 0.075 


' » ° /° 




/^vpc 


-0.044 


± 0.029 


-0.044 


± 0.093 






Ad 


-0.225 


± 0.127 


-0.316 


± 0.199 







ftp 


0.019 


± 0.009 


0.016 


± 0.017 



10 12 10 13 
M vlr (M e ) 

Figure 5. The 3D half-mass radii for stellar masses at the bary- 
onic object scale versus halo mass for the ELO sample. Symbols 
are as in previous Figures 



10" 

Mfir (m ) 



Figure 6. The c rp 



- I fjBtav ra ^; os versus the stellar masses 

3' e,bo 



at the baryonic object scale. Symbols are as in previous Figures. 



grows with decreasing ELO mass scale, and no appreciable 
SF effect. 

Concerning sizes, in Figure 9 we plot the c r d = 
r e°h/ r l*bo ratios versus the Mjf ar mass scale for ELOs in 
both the SF-A and the SF-B samples. In this Figure the ef- 
fects of SF parameterisation are clear: SF-A type ELOs have 
larger sizes relative to the halo size than SF-B type ELOs. 
There is also a clear mass effect, with more massive ELOs 
less concentrated relative to the total mass distribution than 
less massive ones (i.e., spatial homology breaking; note, how- 



Column 2: the slopes of the = MiK® + Mq relation (direct 
fits); the slopes of the Af vir /Afg* ar and c ; tx (M^ T )& scaling 
relations for the the SF-A sample, calculated in log — log plots 
through direct fits. Column 3: their respective 95% confidence 
intervals. Columns 4 and 5: same as columns 2 and 3 for the 
SF-B sample. 




Figure 7. The Af v i r /M^ ar ratios as a function of the ELO mass 
scale. Symbols are as in previous Figures 



ever that the scatter is important). Moreover, Figure 9 sug- 
gests that this trend does not significantly depend on the SF 
parametrisation. These indications are quantitatively con- 
firmed through a fit to a power law c r d = ^rd(M^o ar )' 3rd (see 
Table 5) and have interesting implications to explain the tilt 
of the observed FP. 



4 KINEMATICS 

4.1 Three-dimensional velocity distributions 

Shapes and mass density profiles (i.e., positions) are related 
to the 3D velocity distributions of relaxed E galaxies through 
the Jeans equation (see Binney & Tremaine 1987). Obser- 
vationally, the informations on such 3D distributions is not 
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Figure 8. The M° h /M^ ratios as a function of the ELO mass 
scale. Symbols are as in previous Figures 
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Figure 10. The <T3D(r) profiles of a typical ELO in the SF-A 
sample. Also shown are their anisotropy profiles /3 an i(r). Long- 
dashed lines: dark matter; short-dashed lines: stars 



10 io 10 n 10 is 

Mr (M Q ) 

Figure 9. The c r( j = r*°^ /r st t T ratios as a function of the ELO 
mass scale. Symbols are as in previous Figures. Spatial homology 
breaking is clear in this Figure 



available for external galaxies, only the line-of-sight velocity 
distributions (LOSVD) can be inferred from their spectra. 
They have been found to be close to gaussian (Binney & 
Tremaine 1987; van der Marel & Franx 1993) , so that simple 
equilibrium models can be expected to adequately describe 
their dynamical state (de Zeeuw & Franx 1991). The com- 
plete six dimensional phase space informations for each of 
the particles sampling the ELOs provided by numerical sim- 
ulations, allow us to calculate the velocity profiles, V c ii(r), 



the 3D profiles for the velocity dispersion, 030 (r)> an d their 
corresponding anisotropy profiles. The anisotropy is defined 



/3a 



2(7 2 ' 



(4) 



where a T and at are the radial and tangential velocity 
dispersions (of = cr 2 . + cr^), relative to the centre of the 
object. These profiles, as well as the LOS velocity V\ os (R) 
and LOS velocity dispersion ai OB (R) profiles, are analysed 
in detail in Onorbe et al., in preparation. Figure 10 illus- 
trates the more outstanding results we have obtained: ELO 
velocity dispersion profiles in three dimensions are slightly 
decreasing for increasing r, both for dark matter and stel- 
lar particles, cr^ k (r) and o-f§ r (r). It has been found that 
(°3D k (r)) 2 ~ (1-4 - 2) (afg r (r)) 2 , as Loewenstein 2000 had 
found on theoretical grounds. This is so because stars are 
formed from gas that had lost energy by cooling. This re- 
sult on kinematical segregation is very interesting because 
it has the implication that the use of stellar kinematics to 
measure the total mass of ellipticals could result into inac- 
curate values. The study of anisotropy has shown that it is 
always positive and almost non- varying with r (recall, how- 
ever, that ELOs are non-rotating and that they have been 
identified as dynamically relaxed objects, so that there are 
not recent mergers in our samples). The stellar component 
generally shows more anisotropy than the dark component, 
maybe coming from the radial motion of the gas particles 
that gave rise to the stars. 



4.2 Global parameters for the velocity 
distribution 

Only for a limited number of ellipticals are the V\ OB (R) or 
a\ os {R) profiles available. Observationally, a useful charac- 
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Figure 11. The correlation between the central l.o.s. velocity 
dispersion and the virial mass for the ELO samples. Symbols are 
as in previous Figures. 
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Figure 12. c vpc = (o'| t ^) 2 /3(o-f^ r ) 2 ratios versus the ELO 
mass scale. Symbols are as in previous Figures 



terisation of the velocity dispersion of an E galaxy is pro- 
vided by its central stellar line-of-sight velocity dispersion, 
fflos.O' It corresponds to the velocity dispersion of the stel- 
lar (as opposed to dark matter or other) component. Due 
to its interest, eri os ,o has deserved an important attention 
in literature and it had been measured for several E galaxy 
samples before the SDSS results (Faber et al. 1987; Djorgov- 
ski & Davis 1987; Dressier et al. 1987; Lucey, Bower & Ellis 
1991; j0rgensen, Franx & Kjoergaard 1993; j0rgensen et al. 
1996; Kelson et al. 1997; Kelson et al. 2000, Bernardi et al. 
2002). In Table 4 the values of of£ Q for the ELO sample are 

listed^. In Figure 11 we show the good correlation between 
of*" q and the virial mass, M v i r . 

Physically, a measure of the average dynamical state of 
stars in the ELO itself is provided by their mean square 
velocity relative to the ELO center of mass, or average 
three-dimensional velocity dispersion crf*^, of which the ob- 
servationally available af^ T parameter is assumed to be a 
fair estimator. To test this point, in Figure 12 we plot the 
c vpc = (o"3*bo) 2 /3(o" 1 s *s r ) 2 ratios versus the ELO mass scale. 
We see that no mass effect is apparent, and this is quantita- 
tively confirmed in Table 5, where the results of a fit of the 
form Cvpc = A vpc (M^" r ) l3vpc are given. We also see that due 
to radial anisotropy, c vpc < 1, with no SF parametrisation 
effect. So, there is not mass bias when using af^o as an esti- 
mator for 03 but some warnings are in order concerning 
anisotropy effects. 

A significant velocity dispersion parameter for ELOs 
at the halo scale is 03°^, the average 3-dimensional velocity 
dispersion of the whole elliptical, including both dark and 

^ Recall that the central LOS velocity dispersion for the stellar 
component of ELOs is written as c^J r , with a "star" superindex 
to distinguish it from that of the other ELO components 



baryonic matter. According with Eq. 2, this is the velocity 
dispersion entering in the virial theorem for the whole ELO 
configuration. To test that this is in fact the case, in Fig- 
ure 13 we plot the cp = GM v i T /(al°l) 2 rl°^ form factors (see 
Eq. 2) as a function of M^" 1 . The lack of any significant 
mass or SF parametrisation effects in this Figure are quan- 
titatively confirmed through a fit to power laws of the form 
cp = A-p(M^ I ) l3F , whose results in Table 5 are consistent 
with cp being independent of the ELO mass scale or SF pa- 
rameter values. Note also that the cp values are as expected 
(Binney & Tremaine 1987). 

Once we have confirmed that when writting the virial 
theorem for an ELO configuration, is the velocity dis- 
persion one has to use, let us remind that one has to be 
careful when using af^ T or erf fj* as estimators for this 
physically meaningful quantity. In Figure 14 we plot the 
c v d = (^"h/^ifbo) 2 ratios, that measure how dissipation and 
concentration affect, on average, to the relative values of the 
dispersion at the halo scale (involving also dark matter) and 
at the baryonic object scale. No mass effects are apparent 
in this Figure, but an average kinematical segregation is 
clear, (see Table 5 for the results of a fit to the expression 
Cvd = ^4vd(M^ ar )' 3vd ). These are important results, which 
could have interesting observational implications. 



5 THE INTRINSIC DYNAMICAL PLANES 
AND THE DYNAMICAL FUNDAMENTAL 
PLANES 

5.1 The dynamical plane relations of ELO samples 

In the last section it has been shown that the mass, size 
and velocity dispersion parameters at the halo scale satisfy 
virial relations. This result is, however, at odds with the 
tilt of the observed FP of ellipticals discussed in §1, that 
involves the L, _Rc ght and ci os ,o observed variables, whose 
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Figure 13. The c F EE GM vir / '(o^) 2 ? - * * form factors (see Eq. 2) 
versus the ELO mass scale. Symbols are as in previous Figures. 
This Figure confirms that r*°* and are the size and veloc- 

ity dispersion ELO parameters that must be used in the virial 
theorem 




10" 

Figure 14. The c v( j EE (cg^/Cg*^) 2 ratios (average kinematical 
segregation) as a function of the ELO mass scale. Symbols are 
as in previous Figures. No dynamically broken homology can be 
seen in this Figure, but the kinematical segregation between dark 
matter and stars is clear. 



virtual counterparts describe the ELO at the scale of the 
baryonic object. So, we have first to analyse whether or not 
the mass, size and velocity dispersion of ELOs at this scale 
define planes tilted relative to the virial one. 

To this end, we have carried out a principal component 
analysis (PCA) of the SF-A and SF-B samples in the three 
dimensional variables E = log 10 M^ 1 , r = log 10 JVbo an d 
v = log 10 o-g'j^, through their 3x3 correlation matrix C. 
We have used three dimensional variables rather than pro- 
jected ones to circumvent projection effects, that add noise 
to the correlations (see, for example, Figure 12). We have 
found that, irrespective of the SF parametrisation, one of the 
eigenvalues of C is considerably smaller than the others (see 
Table 6), so that ELOs populate in any case a flattened el- 
lipsoid close to a two-dimensional plane in the (E, r, v) space 
that we call the intrinsic dynamical plane (IDP); the FP is 
the observed manifestation of this IDP. The eigenvectors of 
C indicate that the projection 

E- E= a 3D (r- r) + 7 3D (u- v), (5) 

where E, f and v are the mean values of the E, r and 
v variables, shows the IDP viewed edge-on. Table 6 gives 
the eigenvalues of the correlation matrix C (Ai, A2, A3), the 
planes eq. (5), as well as their corresponding thicknesses 
o~Erv, both for the SF-A and SF-B samples. The IDPs are 
in fact tilted relative to the virial plane (characterized by 
a = 1,7 = 2), and their scatter is very low as measured by 
their thicknesses o~Erv. Note that the eigenvalues of the PCA 
analysis are not dependend on the SF parametrisation. 

In Figure 15 we plot the (E,v), (E,r) and (r, v) projec- 
tions of the IDPs corresponding both to the SF-A sample 
and the SF-B sample. We see that the three projections 
show correlations and that these are very tight for the first 
of them. 



5.2 Comparing the IDP to the observed FP of 
elliptical galaxies 

The next step is to compare the results on the IDPs we 
have found with the tilt and the scatter of the observed 
FP. We note that the r and v variables are not observation- 
ally available, so that we have to use their projected coun- 
terparts. Assuming that the projected stellar mass density 
profile, E star (i?), can be taken as a measure of the surface 
brightness profile, then < E star > c = c < J light > e , with c a 
constant, and -R^bo — ^e' sht an d we can look for a funda- 
mental plane (hereafter, the dynamical FP) in the 3-space of 
the structural and dynamical parameters -R^'boi < E star > c 
and ofoJo) directly provided by the hydrodynamical simula- 
tions. To make this analysis as clear as possible, we trans- 
form to a «;-like orthogonal coordinate system, the dynami- 
cal system, i=l,2,3, similar to that introduced by Bender, 
Burstein & Faber (1992), but using R s c % instead of ^ ight 
and < E star > e instead of < I hsht > c , and, consequently, 
free of age, metallicity or IMF effects. The dynamical nf 
variables can be written as: 

**? = [l°g (ofoM)) 2 + log Rl%] /V2 (6) 
k% = [log (af ta o) 2 + 2 log<E star } - log R B X]/V6 (7) 

K? = [log (af ta 0) 2 - log<E star > e - log Rl%] /y/3 (8) 

and they are related to the k coordinates through the 
expressions: ki — K] 3 , k.2 — kJ? — v / 6/31og(M^ ar /Z/) and 
k 3 ~ k° + v / 3/31og(Mb t a 7L). We discuss the tilt and the 
scatter of the dynamical FP separately. We first address the 
tilt issue. We use at this stage for the -R^bo an( A ^fos.o vari- 
ables the averages over three orthogonal l.o.s. projections, 
to minimize the scatter in the plots caused by projection 
effects. 
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Table 6. Results of PCA at z=0 



Sample 


No. 


E 


f 


V 


Ai 


A 2 


As 


a 3D 


7 3D 


CErv 


SF-A 
SF-B 


26 
17 


10.993 
11.259 


0.746 
0.695 


2.335 
2.453 


0.12893 
0.09342 


0.00355 
0.00292 


0.00013 
0.00014 


0.427 
0.249 


2.066 
2.449 


0.011 

0.012 



Column 2: ELO number in the sample. 
Columns 3, 4 and 5: sample mean values of the E,r and v variables. 
Columns 6, 7 and 8: eigenvalues of the correlation matrix. 
Columns 9 and 10: coefficients of the plane (Eq. (5)). 
Column 11: IDP scatter in the E,r and v variables. 



b 




Figure 15. The IDPs for the SF-A and SF-B samples. Projections on the (E,v), (E,r) and (r,v) coordinate planes are shown. 
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Figure 16. Dynamical Fundamental Plane in k d system. Edge-on projection (top panel) and nearly face-on projection (bottom panel) 
of the dynamical FP of ELOs in the kP variables (black filled symbols: SF-A sample; red open symbols: SF-B sample). We also draw the 
respective concentration ellipses (with their major and minor axes) for the SDSS early-type galaxy sample from Bernardi et al. (2003b) 
in the z band (solid line) and the r band (dashed line). See text for more details. 



Figure 16 plots the k§ versus (top) and k§ ver- 
sus «i (bottom) diagrams for ELOs in both the SF-A and 
SF-B samples. We also drew the 2a concentration ellipses 
in the respective variables, as well as its major and minor 
axes, for the SDSS early-type galaxy sample in the SDSS z 
band (solid lines) and in the r band (point lines) as anal- 
ysed by Bernardi et al. 2003b, 2003c. The most outstanding 
feature of this Figure (upper panel) is the good scaling be- 
haviour of k§ versus Acf , with a very low scatter (see the 
slopes Mi in Table 5; note that the slopes for the SF-A and 
SF-B samples are consistent within their errors, while the 
zero-points depend on the SF parameterisation through the 
ELO sizes) . Another interesting feature of Figure 16 is that 
it shows that most of the values of the ref coefficients are 
within the 2a concentration ellipses in both plots for ELOs 
formed in SF-A type simulations, with a slightly worse agree- 



ment for ELOs in the SF-B sample. This means that ELOs 
have counterparts in the real world (Saiz et al. 2004). It is 
worth mentioning that these results are stable against slight 
changes in the values of the Qa, Hbaryon or h parameters; for 
example, we have tested that using their preferred WMAP 
values (Spergel et al. 2003) shows results negligibly different 
to those plotted in Figure 16. Finally, we note that either 
the dynamical or the observed FPs are not homogeneously 
populated: both SDSS ellipticals and ELOs occupy only a 
region within these planes (see Figure 16 lower panel, see 
also Guzman et al. 1993; Marquez et al. 2000). This means 
that, from the point of view of their structure and dynamics, 
ELOs are a two-parameter family where the two parameters 
are not fully independent. Moreover, concerning ELOs, the 
occupied region changes when the SF parameters change. 
The reason of this change is that the ELO sizes decrease 
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as SF becomes more difficult, because the amount of dis- 
sipation experienced by the stellar component along ELO 
assembly increases. 

We now turn to consider the scatter of the dynamical 
FP for the ELO samples and compare it with the scatter 
of the FP for the SDSS elliptical sample, calculated as the 
square root of the smallest eigenvalue of the 3x3 covariance 
matrix in the E (or logL), V = logoff and R = log 
variables (Saglia et al. 2001). As Figure 16 suggests, when 
projection effects are circumvented by taking averages over 
different directions, the resulting three dimensional orthog- 
onal scatter for ELOs is smaller than for SDSS ellipticals 
(ctevr = 0.0164 and oevr = 0.0167 for the SF-A and SF-B 
samples, respectively, to be compared with olvr. = 0.0489 
for the SDSS in the log L,V = log <7i os ,o and R = log 7?i ight 
variables). To estimate the contribution of projection effects 
to the observed scatter, we have calculated the orthogonal 
scatter for ELOs when no averages over projection direc- 
tions for the i?o'bo an d o"fos,o variables are made. The scatter 
(vevr = 0.0238 and o E vr = 0.0214 for the SF-A and SF-B 
samples) increases, but it is still lower than observed. This 
indicates that a contribution from stellar population effects 
is needed to explain the scatter of the observed FP, as sug- 
gested by different authors (see, for example, Pahre et al. 
1998; Trujillo et al. 2004). 

5.3 Clues on the physical origin of the IDP tilt 

We now address the issue of the physical origin of the tilt of 
ELO IDPs relative to the virial relation. As discussed in §1, 
a non-zero tilt can be caused by a mass dependence of the 
mass-to- light ratio M v i r /L, of the mass structure coefficients 
Cu = — Gi^vir Qr Q £ koth of them. We examine these 

M 3d-, 2 „«' 

possibilities in turn. 

i) We first note that the mass-to-light ratio can be writ- 
ten as: 

M vil /L = A vil (MiTf-" x 7 star , (9) 

where 7 star = M^ r / L is the stellar mass-to-light ratio, 
that, as already explained, can be considered to be indepen- 
dent of the E galaxy luminosity or ELO mass scale. Figure 7 
and the values of the /3 v i r slopes given in Table 5, indicate 
that the dark to bright mass content of ELOs increases with 
their mass, contributing a tilt to their IDPs. Similar results 
have also been found in pre-prepared simulations of dissipa- 
tive mergers (Robertson et al. 2006). 

ii) Writting the c M r mass structure coefficients as power 
laws c M r = AMiM^ 1 ) 13 ^, ELO homology would imply /3 M = 
0. To elucidate whether or not this is the case, the /3m slopes 
have been measured on the ELO samples through direct fits 
in log-log scales. The results are given in Table 5, where we 
see that the homology is in fact broken both for SF-A or 
SF-B samples. To deepen into the causes of this behaviour, 
we use Eqs. 2 and 3 to write 

C M — CpC r dC rp C v dCvpc, (10) 

where the c\ coefficients, with i = rd, rp, vd, vpc have 
been defined in §3.2, and with i = F in §4.2. Taking into 
account the power-law forms of these coefficients, we have: 

/3m =/3F+/3 r d+/3rp+/3vd+/3vpc, (11) 



when the /3i slopes are calculated through direct fits. 
From Table 5 we see that the main contribution to the ho- 
mology breaking comes from the c r d coefficients (i.e., spatial 
homology breaking, see §3.2 and Onorbe et al. 2005), while, 
as previously discussed, no dynamical homology breaking 
(see §4.2) or projection effects are important in our ELO 
samples. 



6 SUMMARY, DISCUSSION AND 
CONCLUSIONS 

6.1 Summary 

We present an analysis of the sample of ELOs formed in 
ten different cosmological simulations, run within the same 
global flat A cosmological model, roughly consistent with 
observations. The normalisation parameter has been taken 
slightly high, as = 1.18, as compared with the average fluc- 
tuations of 2dFGRS or SDSS galaxies to mimic an active 
region of the Universe. Newton laws and hydrodynamical 
equations have been integrated in this context, with a stan- 
dard cooling algorithm and a star formation parameteriza- 
tion through a Kennicutt-Schmidt-like law, containing our 
ignorance about its details at sub-kpc scales. No further hy- 
potheses to model the assembly processes have been made. 
Individual galaxy-like objects naturally appear as an output 
of the simulations, so that the physical processes underlying 
mass assembly can be studied. Five out of the ten simula- 
tions (the SF-A type simulations) share the SF parameters 
and differ in the seed used to build up the initial conditions. 
To test the role of SF parameterisation, the same initial con- 
ditions have been run with different SF parameters making 
SF more difficult, contributing another set of five simulations 
(the SF-B type simulations). ELOs have been identified in 
the simulations as those galaxy-like objects at z — having 
a prominent, non-rotating dynamically relaxed spheroidal 
component made out of stars, with no extended discs and 
very low gas content. These stellar component is embedded 
in a dark matter halo that contributes an important fraction 
of the mass at distances from the ELO centre larger than 
~ 10-15 kpc. No ELOs with stellar masses below 3.8 x 10 10 
Mq or virial masses below 3.7 x 10 11 M have been found 
that met the selection criteria (see KaufTmann et al. 2003b 
for a similar result in SDSS galaxies and Dekel & Birnboim 
2006, and Cattaneo et al. 2006 for a possible physical expla- 
nation). ELOs have also an extended halo of hot, diffuse gas. 
Stellar and dark matter particles constitute a dynamically 
hot component with an important velocity dispersion, and, 
except in the very central regions, a positive anisotropy. 

The informations about position and velocity distribu- 
tions of the ELO particles of different kinds (dark matter, 
stars, cold gas, hot gas) provided by the simulations, allows 
a detailed study of their intrinsical three dimensional mass 
and velocity distributions, as well as a measure of the pa- 
rameters characterising their structure and dynamics. In a 
forthcoming paper, we report on the three dimensional mass 
density, circular velocity and velocity dispersion profiles, as 
well as the projected stellar mass density profiles and the 
LOS velocity dispersion profiles. In this paper we focus on a 
parameter analysis to quantify some of the previous results. 

Mass, size and velocity dispersion scales for their differ- 
ent components have been measured in the ELO samples, 
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both at the scale of their halo and at the scale of the bary- 
onic object (a few tens of kiloparsecs) . At the halo scale, the 
masses of both cold gas and stars, M^ h and M{f ar , respec- 
tively, have been found to be tightly correlated with the halo 
total mass, M vir , with the ratios M^ h /M vir and M^ al /M vil 
decreasing as A/ vir increases (that is, massive objects miss 
cold baryons within r v j r when compared with less massive 
ELOs) , presumably because gas gets more difficulties to cool 
and fall as M vir increases. The overall half-mass radii, r^h 
shows also a very tight correlation with M v ir • Half-mass radii 
for the cold baryon or stellar mass distributions have a more 
complex behaviour, as in these cases gas heating in shocks 
and energy losses due gas cooling are in competition to de- 
termine these distributions. 

A very important result we have found when analysing 
ELOs at the scale of the baryonic object, is that M v i r plays 
an important role to determine the ELO structure also 
below a few tens of kiloparsecs scales. In fact, both the 
masses of cold baryons M^ (i.e., those baryons that have 
reached the central regions of the configuration), and of stars 



M£ 



show a good correlation with M v i r , and, moreover, 



the M^/Mvir and M^ 1 /M vil ratios (i.e., the relative con- 
tent of cold baryons or stars versus total mass) decrease 
as M v i r increases. This is the same qualitative behaviour 
shown by these ratios observationally in the SDSS data, 
and, also, by ELOs at the halo scale. The dependence of 

or M^o" on the SF parametrisation is only very slight, 
with SF-A type ELOs having slightly more stars than their 
SF-B type counterparts. The half-mass radii for cold baryon 
and stellar masses, rg h ho and r^ T a , show also a good cor- 
relation with M v i r , but now the values of the SF parame- 
ters also play a role, because their change implies a change 
in the time interval during which gas cooling is turned on, 
and this changes the ELO stellar mass distribution, i.e., its 
lenghthscale, so that ELO compactness increases from SF-A 
to SF-B type simulations. Another important result is that, 
regardless of the SF parametrizations used in this work, the 
relative distributions of the stellar and dark mass compo- 
nents in ELOs show a systematic trend measured through 
the Crd = fc°h/ r "c*bo ratios, with stars relatively more concen- 
trated as A/ V ir decreases (i.e., a quantification of the spatial 
homology breaking). Note that to compare with observa- 
tional data, the relevant parameters are the projected half- 
mass radii, -R^w We have checked that they show an ex- 
cellent correlation with the corresponding three dimensional 
half-mass radii, with the c rp = r^ /Rf^ ratios showing no 
significant dependence on the ELO mass scale. 

Concerning kinematics, a useful characterisation of the 
ELO velocity dispersion is the central stellar line-of-sight 
velocity dispersion, of g r , whose observational counterpart 
can be measured from elliptical spectra. A very important 
result is the very tight correlation we have found between 
Mvir and ofoJo, confirming that the observationally measur- 
able <ti OSj o is a fair virial mass estimator. In addition, af^* 
is closely related to the mean square velocity of both, the 
whole elliptical at the halo scale (including the dark mat- 



ter), 



'3, hi 



and the stellar component of the central object, 



star 



'3, bo- We have also found that the c v d = ( CT 3^h/ CT 3*bo) or the 
c vp c = (o"3*bo) 2 /3(cr 1 s *" ) 2 ratios are roughly independent of 
the ELO mass scale. And so, ELOs do not show dynamically 
broken homology, even if their stellar and dark components 
are kinematically segregated (i.e., c v d 7^ 1). This could lead 



M£r(M ) 



Figure 17. The M^/M* ratios as a function of the ELO mass 
scale for three ELOs identified at Z = 1 in S128 (red) and their 
counterparts in S64 (black). 



to inaccurate determinations of the total mass of ellipticals 
when using stellar kinematics. 

A very important result is that, irrespective of the SF 
parameterisation, the (logarithms of the) ELO stellar masses 
A/si* ar , g^gjjgj. half-mass radii r^ T a , and stellar mean square 
velocity of the central object 03*^, define intrinsic dynam- 
ical planes (IDPs). These planes are tilted relative to the 
virial plane and the tilt does not significantly depend on 
the SF parameterisation, but the zero point does depend. 
Otherwise, the intrinsic dynamical plane is not homoge- 
neously populated, but ELOs, as well as E galaxies in the 
FP (Guzman et al. 1993), occupy only a particular region 
defined by the range of their masses. 



6.2 Testing possible resolution effects 

To make sure that the results we report in this paper are 
not unstable under resolution changes, a control simulation 
with 128 3 dark matter and 128 3 baryonic particles, a gravi- 
tational softening of e = 1.15 kpc and the other parameters 
as in SF-A type simulations (the S128 simulation), has been 
run. The results of its analysis have been compared with 
those of a 2 x64 3 simulation (the S64 simulation), whose 
initial conditions have been built up by randomly choosing 
1 out of 8 particles in the S128 initial conditions, so that 
every ELO in S128 has a counterpart in the lower resolution 
simulation and conversely. Due to the very high CPU time 
requirements for S128, the comparison has been made at 
z = 1. The results of this comparison are very satisfactory, 
as Figure 17 illustrates (compare with Figure 7). 
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6.3 The dimensionality of ELO and elliptical 
samples in parameter space 

The intrinsic dynamical planes and their occupations (see 
§ 6.1) reflect the fact that dark matter haloes are a two- 
parameter family (for example, the virial mass and the en- 
ergy content or the concentration; see, for example, Hern- 
quist 1990; Navarro, Frenk & White 1995, 1996; Manrique et 
al. 2003; Navarro et al. 2004) where the two parameters are 
correlated (see, for example, Bullock et al. 2001; Wechsler 
et al. 2002; Manrique et al. 2003). Adding gas implies that 
heating and cooling processes also play a role at determin- 
ing the mass and velocity distributions, and, more particu- 
larly, the length scales. However, as explained above, we have 
found that both, the relative content and the relative dis- 
tributions of the dark and baryonic mass components show 
systematic trends with the ELO mass scale, that can be writ- 
ten as power-laws of the form M vir /M^ ar = A vir (Mf ar ) /3 ™ 
and ri%/rl% = A vd (MZr) f3 < d - 

A first consequence of the regularity of the trends with 
the mass scale found in this paper is that no new parameters 
are added relative to the dark matter halo family, so that the 
baryonic objects are also a two-parameter family, and ELO 
structural and dynamical parameters define also a plane. A 
second consequence is that the plane is tilted relative to the 
halo plane (i.e., the virial plane) because /3 r< j — f3 v i T 7^ 0. 
Finally, the plane is not homogeneously populated because 
of the mass-concentration halo correlation, that at the scale 
of the baryonic objects appears for example as a mass — size 
correlation. This explains the role played by M V [ T to de- 
termine the intrinsic three dimensional correlations. In this 
paper we also show that of^^o is a fair empirical estimator 
of Mvir, and this explains the central role played by cri os ,o at 
determining the observational correlations. 

The fundamental plane shown by real elliptical samples 
is the observationally manifestation of the IDPs when using 
projected parameters Rf^, crf^ and luminosity variables 
instead of stellar masses M^f 1 • We have taken advantage of 
the constancy of the stellar-mass-to-light ratios of ellipticals 
in the SDSS (Kauffmann et al. 2003a, 2003b) to put the 
elliptical sample of Bernardi et al. (2003b, 2003c) in the same 
projected variables we can measure in our virtual ellipticals. 
We have found that the FPs shown by the two ELO samples 
are consistent with that shown by the SDSS elliptical sample 
in the same variables, with no further need for any relevant 
contribution from stellar population effects to explain the 
observed tilt. These effects could, however, have contributed 
to the scatter of the observed FP, as the IDPs have been 
found to be thinner than the observed FP. 



6.4 The physical origin of the tilt in a 
cosmological context 

We now turn to discuss the physical origin of the trends 
given by the power laws M vir /Mb* ar = A^M^) 13 "'" and 
<h/<bo = A lA {M^f^. As explained in Section 2, the 
simulations provide us with clues on the physical processes 
involved in elliptical formation (see also DSS04 and DTal06). 
They also indicate that the dynamical plane appears at an 
early violent phase as a consequence of ELO assembly out 
of gaseous material, with cooling and on short timescales, 
and it is preserved during a later, slower phase, where dis- 



sipationless merging plays an important role in stellar mass 
assembly (see more details in DTal06) . Our simulations show 
that the physical origin of the trends above lie in the sys- 
tematic decrease, with increasing ELO mass, of the rel- 
ative amount of dissipation experienced by the baryonic 
mass component along ELO stellar mass assembly (DTal06; 
Onorbe et al., in preparation). This possibility had been sug- 
gested by Bender et al. (1992), Guzman et al. (1993) and 
Ciotti et al. (1996). Bekki (1998) first addressed it numeri- 
cally in the framework of the merger hypothesis for ellipti- 
cal formation through pre-prepared simulations of dissipa- 
tive mergers of disk galaxies, where the rapidity of the star 
formation in mergers is controlled by a free efficiency param- 
eter Csf- He shows that the star formation rate history of 
galaxies determine the differences in dissipative dynamics, 
so that to explain the slope of the FP he needs to assume 
that more luminous galaxies are formed by galaxy mergers 
with a shorter timescale for gas transformation into stars. 
Recently, Robertson et al. (2006) have confirmed these find- 
ings on the importance of dissipation to explain the FP tilt. 

In this paper we go an step further and analyze the 
FP of virtual ellipticals formed in a cosmological context, 
where individual galaxy-like objects naturally appear as an 
output of the simulations. Our results essentially include 
previous ones and add important new informations. First, 
our results on the role of dissipation to produce the tilt of 
the FP essentially agree with those obtained through dis- 
sipative pre-prepared mergers, but it is important to note 
that, moreover, more massive objects produced in the sim- 
ulations do have older means and narrower spreads in their 
stellar age distributions than less massive ones (see details 
DSS04); this naturally appears in the simulations and need 
not be considered as an additional assumption. Second, the 
preservation of the FP in the slow phase of mass aggre- 
gation in our simulations also agrees with previous work 
based on dissipationless simulations of pre-prepared merg- 
ers (Capelato et al. 1995; Dantas et al. 2003; Gozalez-Garcfa 
& van Albada 2003; Boylan-Kolchin et al. 2005; Nipoti, Lon- 
drillo & Ciotti 2003) . Moreover, elliptical properties recently 
inferred from observations (for example, the appearence of 
blue cores, Menanteau et al. 2004, and the increase of the 
stellar mass contributed by the elliptical population since 
higher z, Bell et al. 2004; Conselice, Blackburne, & Papovich 
2005; Faber et al. 2005; see more details in DTal06) can also 
be explained in our simulations. 

6.5 Conclusions 

We conclude that the simulations provide an unified scenario 
where most current observations on ellipticals can be inter- 
related. In particular, it explains the most important results 
relative to the physical origin of their FP relation (i.e. the FP 
tilt is due dissipative dynamics, and disipationless merging 
in the slow growth phase of mass assembly does not change 
the tilt). It is worth mentioning that this scenario shares 
some characteristics of previously proposed scenarios, but 
it has also significant differences, mainly that most stars in 
elliptical galaxies form out of cold gas that had never been 
shock heated at the halo virial temperature and then formed 
a disk, as the conventional recipe for galaxy formation pro- 
pounds (see discussion in Keres et al. 2005 and references 
therein). The scenario for elliptical formation emerging from 
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our simulations has the advantage that it results from simple 
physical laws acting on initial conditions that are realiza- 
tions of power spectra consistent with observations of CMB 
anisotropies. 
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